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X-RAY EMISSION FROM THE WOLF-RAYET BUBBLE S 308* 
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ABSTRACT 

The Wolf-Rayet (WR) bubble S 308 around the WR star HD 50896 is one of the only two WR bubbles 
known to possess X-ray emission. We present XMM-Newton observations of three fields of this WR 
bubble that, in conjunction with an existing observation of its Northwest quadrant, map most of the 
nebula. The X-ray emission from S 308 displays a limb-brightened morphology, with a central cavity 
~22' in size and a shell thickness of ~8'. This X-ray shell is confined by the optical shell of ionized 
material. The spectrum is dominated by the He-like triplets of N VI at 0.43 keV and O VII at 0.57 keV, 
and declines towards high energies, with a faint tail up to 1 keV. This spectrum can be described by a 
two-temperature optically thin plasma emission model (Ti ~1.1 xlO 6 K, T2 -~13xl0 6 K), with a total 
X-ray luminosity ~ 2xl0 33 erg s _1 at the assumed distance of 1.5 kpc. 

Subject headings: ISM: bubbles - ISM: individual (S308) - stars: individual (HD 50896) - stars: 
winds, outflows - stars: Wolf-Rayet - X-rays: individuals 308) 



1. INTRODUCTION 

Wolf-Rayet (WR) bubbles are the final result of the 
evolution of the circumstellar medium (CSM) of massive 
stars with initial masses M >35 Mq. These stars ex- 
hibit high mass-loss rates throughout their lives, peaking 
during their post-main-sequence evolution that involves 
a Red or Yellow Supergiant (RSG or YSG) or Lumi- 
nous Blue Variable (LBV) stage (e.g.. lMevnet fc Maederl 
12003ft during which the mass-loss rate can be as high 
as 10- 4 -10" 3 M Q yr" 1 (jNugis fc Lamersl2000ft . although 
the stellar wind velocity is low (10-10 2 km s _1 ). The final 
WR stage is characterized by a fast stellar wind (voo > 
10 3 km s _1 ), which sweeps up, shocks, and compresses 
the RSG/LBV material. Thin-shell and Rayleigh- Taylor 
instabilities lead to the corrugatio n and eventual frag- 
menta ti on of the swept-up she l l (iGarcia-Segura et alj 
[l996alrbl : IFrever et al .1 l200l 120061: IToala fc Arthurll2011h . 
Clumpy WR wind-blown bubbles have been detected at 
optical wavelengths around ~10 WR stars in our Galaxy 
(IChu et al.l H983I : IGruendl et all [20001 IStock fc Barlow! 
12010ft . Their optical emission is satisfactorily mod- 
eled as photoionized dense clumps and shell material 
(jEsteban et al.lll993ft . 

X-ray emission has been detected so f ar in only 
two WR bubbles NGC 6888 and S308 (IBochkarev 
1988 1 IWrigge et al.lll994l. [j~998l: IWriggdll99fl: IChu et al 



20031 iVvrigge fc Wendkerl 120021: IWrigge" et a. . I 120051 : 
Zhekov fc Park! 12011ft . The most sensitive X-ray ob- 



servations of a WR bubble are thos e of the no r thwes t 
(NW) quadrant of S 308 presented by IChu et ail (|2003ft . 
Their XMM-Newton EPIC-pn X-ray spectrum of S 308 
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revealed very soft X-ray emission dominated by the N VI 
He-like triplet at -~0.43 keV and declining sharply to- 
ward higher energies. This spectrum was fit with a 
two-temperature optically thin MEKAL plasma emission 
model, with a cold main component at fcTi = 0.094 keV 
(i.e., Tx ~ 1.1 x 10 6 K), and a hot secondary component 
at kT-i ~ 0.7 keV contributing <6% of the observed X-ray 
flux. The comparison of the X-ray and optical Ha and 
[O III] images of S 308 showed that the X-ray emission is 
confined by the ionized shell. 

In this paper, we present the analysis of three ad- 
ditional XMM-Newton observations of S308, which, in 
co njunction with th ose of the NW quadrant presented 
by IChu et all (|2u03ft . map 90% of this WR bubble (see 
§2). In §3 and §4 we discuss the spatial distribution and 
spectral properties of the X-ray-emitting plasma in S 308, 
respectively. In §5 we present our results of the X-ray 
emission from the central star in the WR bubble. A dis- 
cussion is presented in §6 and summary and conclusions 
in §7. 

2. XMM-NEWTON OBSERVATIONS 

The unrivaled sensitivity of the XMM-Newton EPIC 
cameras to large-scale diffuse emission make s them the 
prefer red choice for the observation of S 308. IChu et all 
(|2003ft presented XMM-Newton observation s of the 
brightest NW quadrant of S 308 (|Wriggel 1 1999ft . but the 
large angular size of S308 (~40/ in diameter) exceeds 
the field of view of the EPIC camera and a significant 
fraction of the nebula remained unobserved. To comple- 
ment these observations, additional XMM-Newton obser- 
vations of three overlapping fields covering the northeast 
(NE), southwest (SW), and southeast (SE) quadrants of 
the nebula have been obtained. The previous and new 
observations result in the coverage of ~90% of the area 
of S 308. The pointings, dates and revolutions of the ob- 
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Fig. 1. — XMM-Newton EPIC images of the four observations of S308 in the 0.3—1.15 keV band. The images have been extracted using 
a pixel size of 2'.'0 and adaptively smoothed using a Gaussian kernel between 5" and 30". The source regions used for spectral analysis 
are indicated by solid lines and the background regions by dotted lines. Note that the point sources that are present in the images were 
excised for the spectral analysis. 



servations, and their exposure times are listed in Table[TJ 
In the following, we will refer to the individual observa- 
tions by the quadrant of S 308 that is covered, namely 
the NW, NE, SW, and SE quadrants. All observations 
were obtained using the Medium Filter and the Extended 
Full-Frame Mode for EPIC-pn and Full-Frame Mode for 
EPIC-MOS. 

The XMM-Newton pipeline products were processed 
using the XMM-Newton Science Analysis Software (SAS) 
Version 11.0, and the Calibration Access Layer avail- 
able on 2011-09-13. In order to analyze the diffuse and 
soft X-ray emission from S308, the XMM-Newton Ex- 
tended Source Analysis Software (XMM-ESAS ) package 
(jSnowden et al.ll2008b iKuntz fc Snowdenl[2008h has been 



used. This procedure applies very restrictive selection 
criteria for the screening of bad events registered during 
periods of high background to ensure a reliable removal 
of the background and instrumental contributions, par- 
ticularly in the softest energy bands. As a result, the 
final net exposure times resulting from the use of the 
XMM-ESAS tasks, as listed in Table [TJ are noticeably 
shorter than the original exposure times. Since we are 
interested in the best time coverage of the central WR 
star to assess its possible X-ray variability and given that 
its X-ray emission level is much brighter than that of 
a mildly enhanced background, we applied less restric- 
tive criteria in selecting good time intervals for this star. 
For this particular analysis, the 10-12 keV energy band 
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TABLE 1 

XMM-Newton Observations of S 308 



Pointing 


R.A. Dec. 


Obs. ID 


Rev. 


Observation start 


Total 


exposure 


time 


Net 


exposure 


time 




(J2000) 






UTC 


pn 


MOS1 


MOS2 


pn 


MOS1 


MOS2 












[ks] 


[ks] 


[ks] 


[ks] 


[ks] 


[ks] 


NW 


06:53:30 -23:43:00 


0079570201 


343 


2001-10-23T22:00:09 


43.5 


47.6 


47.5 


11.9 


19.6 


19.9 


SW 


06:53:24 -23:56:18 


0204850401 


781 


2004-03-15T14:30:41 


20.0 


23.3 


23.4 


6.4 


9.0 


9.2 


SE 


06:55:16 -24:00:00 


0204850501 


781 


2004-03-14T23:00:41 


22.0 


25.4 


25.4 


8.2 


12.4 


12.7 


NE 


06:54:47 -23:46:18 


0204850601 


781 


2004-03-15T06:45:41 


22.0 


25.4 


25.4 


5.4 


8.9 


8.4 



is used to assess the charged particle background, and 
we excised periods of high background with EPIC-pn 
count rates >1.5 counts s~ : and EPIC-MOS count rates 
>0.3 counts s" 1 . 

3. SPATIAL DISTRIBUTION OF THE DIFFUSE 
X-RAY EMISSION 

3.1. Image processing 

Following Snowden & Kuntz' cookbook for anal- 
ysis procedures for XMM-Newton EPIC obser- 
vations of extended ob jects and diffus e back- 
ground, Version 4.3 (Snowden fc Kuntzl |201~TI ). 
the XMM-ESAS tasks and the associated Cur- 
rent Calibration Files (CCF), as obtained from 
|f tp: //xinni. esac . esa. int/pub/ ccf /const ituents/extras/esas. 
have been used to remove the contributions from the 
astrophysical background, soft proton background, 
and solar wind charge-exchange reactions, which have 
contributions at low energies (<1.5 keV). The resulting 
exposure-map-corrected, background-subtracted EPIC 
images of each observed quadrant of S 308 in the 
0.3-1.15 keV band are presented in Figure [TJ The new 
observations of the NE, SW, and SE quadrants of S 308 
detect diffuse emission, as well as a significant number 
of point sources superimposed on this diffuse emission. 
With the single exception of HD 50896 (a.k.a. WR6), 
the WR star progenitor of this bubble registered in the 
SW and NE observations, all point sources are either 
background or foreground sources that we have removed 
prior to our analysis. 

3.2. Analysis of the diffuse X-ray emission 

In order to analyze the spatial distribution of the dif- 
fuse X-ray emission in S 308, the four individual observa- 
tions have been mosaicked using the XMM-ESAS tasks 
and all point sources removed using the Chandra Interac- 
tive Analysis of Observations (CIAO) Version 4.3 dmfilth 
routine, except the one corresponding to WR6. The fi- 
nal image (Figure HJZe/f), extracted in the 0.3-1.15 keV 
energy band with a pixel size of 3'.'0, has been adaptively 
smoothed using the ESAS task adapt-2000 requesting 
100 counts of the original image for each smoothed pixel, 
with typical smoothing kernel scales ^ 1' in the bright- 
est regions and 1' — 2' in the faintest ones. This image 
reveals that the diffuse X-ray emission from S 308 has a 
limb-brightened morphology, with an irregular inner cav- 
ity ~ 22' in size. The surface brightness distribution dis- 
played by this image confirms and adds fur ther details to 
the results of pr evious X-ray observations (|Wriggelll999l : 
IChu et al.ll2003h . The X -ray emission from the bubble 
is brighter towards the northwest blowout and the west- 
ern rim, and fainter towards the east. The bubble seems 



to lack detectable X-ray emission towards the central re- 
gions around the WR star. 

The limb-brightened spatial distribution of the X-ray 
emission from S 308 is further illustrated by the surface 
brightness profiles along the SE-NW and NE-SW direc- 
tions shown in Figure [3J The emission in the innermost 
regions, close to the central WR star, falls to levels com- 
parable to those of a background region to the west of 
S 308 shown in the bottom panel of Fig. [3] Besides the SE 
region, whose surface brightness distribution is best de- 
scribed by a plateau, the X-ray emission along the other 
directions increases steadily with radial distance, peak- 
ing near the shell rim and declining sharply outwards. 
The thickness of the X-ray-emitting shell is difficult to 
qua ntify: along the SW direction, it has a FWHM ~5', 
- c w}$reas it has a FWHM ~8' along the NE direction. 
Figure [3] also illustrates that the X-ray-emitting shell is 
larger along the SE-NW direction (~44' in size) than 
along the NE-SW direction (~40' in size). 

Finally, the spatial distribution of the diffuse X-ray 
emission from S 308 is compared to the [O III] emission 
from the ionized optical shell in Figure fright. The X- 
ray emission is interior to the optical emission not only 
for the NW quadrant but for the entire bubble. This is 
also illustrated in the color composite picture shown in 
Figure 31 in which the distribution of the X-ray emission 
is compared to the optical Ha and [O III] images. This 
image shows that the diffuse X-ray emission is closely 
confined by the filamentary emission in the Ha line, 
whereas the smooth emission in the [O III] line extends 
beyond both the Ha and X-ray rims. 



4. PHYSICAL PROPERTIES OF THE HOT GAS IN 

S308 

The spectral properties of the diffuse X-ray emission 
from S 308 can be used to investigate the physical 
conditions and chemical abundances of the hot gas 
inside this nebula. In order to proceed with this 
analysis, we have defined several polygonal aperture 
regions, as shown in Figure [TJ which correspond to 
distinct morphological features of S308: regions with 
#1 designations correspond to the rim revealed by the 
limb-brightened morphology, #2 the NW blowout, #3, 
#5, and #6 shell interior, and #4 the central star 
HD 50896. We note that any particular morphological 
feature may have been registered in more than one 
quadrant, in which case several spectra can have the 
same numerical designation (for instance, there are four 
spectra for the rim of the shell, namely 1NE, 1NW, 1SE, 
and 1SW). 
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Fig. 2. — (left) Adaptively smoothed XMM-Newton EPIC Image of S308 in the 0.3-1.15 keV band. All point sources, except for the 
central star HD 50896 (WR6), have been excised, (right) Ground- based [O III] image of S 308 obtained with the Michigan Curtis Schmidt 
telescope at Cerro Tololo Inter-American Observatory (CTIO) with superimposed X-ray emission contours. The position of the central star 
HD 50896 is indicated in both panels. 



4.1. Spectra Extraction and Background Subtraction 

Perhaps the most challenging problem associated with 
the analysis of the X-ray spectra of S 308 is a reliable 
subtraction of the background contribution. The diffuse 
X-ray emission from S 308 fills a significant fraction of 
the field of view of the EPIC-pn and EPIC-MOS cam- 
eras, making the selection of suitable background regions 
difficult because the instrumental spectral response of the 
cameras close to their edges may not be the same as those 
for the source apertures. 

The background contribution to the diffuse emission 
from clusters of galaxies that fills the field of view is 
typically assessed from high signal-to-noise ratio obser- 
vations of blank fields. In the case of S308, however, 
the comparison of spectra extracted from background 
regions with those extracted from the same detector re- 
gions of the most suit able EPIC Blank Sky observations 
(|Carter fc~ Rcad 20Q3) clearly indicates that they have 
different spectral shapes. The reason for this discrep- 
ancy lies in the typical high Galactic latitude of the EPIC 
Blank Sky observations, implying low hydrogen absorp- 
tion column densities and Galactic background emission, 
while S 308 is located in regions close to the Galactic 
Plane where extinction and background emission are sig- 
nificant. We conclude that EPIC Blank Sky observa- 
tions, while suitable for the analysis of the diffuse emis- 
sion of a large variety of extragalactic objects, cannot be 
used in our analysis of S 308. 

Alternatively, the different contributions to the com- 
plex background emission in XMM-Newton EPIC obser- 
vations can be modeled, taking into account the contri- 
butions from the astrophysical background, solar wind 
charge-exchange reactions, high-energy (soft protons) 
particle contributions, and electronic noise. This is the 



procedure recommended by the XMM-ESAS in the re- 
lease of SAS vll.O, followi ng the backgro u nd m odel- 
ing methodology devised by iSnowden et al.l (|2004f ) and 
iKuntz fc Snow dcn (2008[). Even though the modeling of 
the different contributions is a complex task, it can be 
routinely carried out. Unfortunately, S308 is projected 
close to the Galactic Plane and the ROSAT All Sky Sur- 
vey (RASS) reveals that it is located in a region of strong 
soft background emission with small-scale spatial varia- 
tions. As shown in Figure [5J the X-ray emission from 
this background is soft and shows lines in the 0.3-1.0 
keV energy band from thermal components as the emis- 
sion from S 308. Therefore, it is not possible to model 
independently the emission from the S 308 bubble and 
that from the soft background. 

The only viable procedure for the analysis seems the 
use of background spectra extracted from areas near the 
camera edges of the same observations. It can be ex- 
pected that the spectral properties of these background 
spectra differ from those of the background registered 
by the central regions of the cameras, given the varying 
spectral responses of the peripheral and central regions of 
the cameras to the various background components. To 
assess these differences, we have used EPIC Blank Sky 
observations to extract spectra from source and back- 
ground regions identical in detector units to those used 
for S 308. Two typical examples of blank sky background- 
subtracted spectra are presented in Figure|Hl While these 
spectra are expected to be flat, several deviations can be 
noticed: (1) clear residuals at ~ 7.5 and ~ 8.1 keV, which 
can be attributed to the defective removal of the strong 
instrumental Cu lines th at affect the EPIC-pn spectra 
(|Kuntz fc Snowdenl 120081 ). (2) a noticeable deviation at 
~ 0.65 keV of the O-K line, and (3) most notably, de- 
viations at energies below 0.3 keV, which is indicative 
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Fig. 3. — S308 X-ray surface brightness profiles along the SE- 
NW (PA=135°) and SW-NE (PA=45°) directions extracted from 
the smoothed XMM-Newton EPIC image presented in Figure [2} 
left. For comparison, a surface brightness profile of a representative 
background region towards the West of S 308 is shown at the same 
intensity and spatial scales. 

of the faulty remo val of the electronic n oise component 
of the background ([Carter fe Readll2007t) . Note that the 
Al-K line at ~ 1.5 keV and the Si-K line at - 1.8 keV, 
which may be expected to be strong in the background 
EPIC-pn spectra, are correctly removed. Consequently, 
we have chosen to use the background spectra extracted 
from the observations of S308, but restrict our spectral 
fits to the 0.3-1.3 keV band. 

4.2. Spectral properties 

The individual background-subtracted EPIC-pn spec- 
tra of the diffuse emission of S 308 are presented in Fig- 
ure [7J This figure also includes spectra of the whole neb- 
ula, its rim, and central cavity obtained by combining all 
spectra, those of regions #1 and #2 (the rim or limb- 
brightened shell), and those of regions #3, #5, and #6 
(the central cavity), respectively. Spectra and response 
and ancillary calibrations matrices from different obser- 
vations of the same spatial regions were merged using 
the mathpha, addarf, and addrmf HEASOFT tasks ac- 
cording to the prescriptions of the SAS threads. The 
EPIC-MOS spectra have also been examined and found 
to be consistent with the EPIC-pn ones although, due 
to the lower sensitivity of the EPIC-MOS cameras, they 
have fewer counts. Therefore, our spectral analysis of the 
diffuse X-ray emission will concentrate on the EPIC-pn 
spectra. 

The spectra shown in Figure [7] are all soft, with a 



prominent peak near > 0.4 keV, and a rapid decline in 
emission towards higher energies. Some spectra (e.g., 
1NE, 1NW, 1SW, 5NE, and 5SE) show a secondary peak 
near < 0.6 keV that is only hinted in other spectra. There 
is little emission above ~ 0.7 keV. although some spectra 
(e.g., 1SE, 1SW, 2NW, 3NE, 5SE, and 6SW) appear to 
present a hard component between 0.8 and 1.0 keV. 

The feature at ~ 0.4 keV can be identified with the 
0.43 keV N VI triplet, while the fainter feature at ~ 
0.6 keV can be associated with the 0.57 keV O VII triplet. 
The occurrence of spectral lines is suggestive of optically 
thin plasma emissi on, confirming previous X-ray s pectral 
analyses of S308 (|Wriggd[l999[ IChu et al.|[2003l) . The 
predominance of emission from the He-like species of ni- 
trogen and oxygen over their corresponding H-like species 
implies a moderate ionization stage of the plasma. Fur- 
thermore, the relative intensity of the N VI and O VII 
lines suggests nitrogen enrichment, since the intensity of 
the O VII lines from a plasma with solar abundances 
would be brighter than that of the N VI lines. 

In accordance with their spectral properties and 
previous sp ectral fits of the NW regions of S 308 
(|Chu et al.l 120031 ). all the X-ray sp ectra of S308 
have been fit with XSPEC vl2.7.0 (|Arnaudl [1991 
for an absorbed two-temperature APEC optically thin 
plasma emission model. The absorpti on model uses 
iBalucinska-Church fe McCammonl (|1992| ) cross-sections. 
A low temperature component is used to model the 
bulk of the X-ray emission, while a high temperature 
component is used to model the faint emission above 
0. 7 keV. We have ad opted the same chemical abundances 
as IChu et all ([2001 ), i.e., C, N, O, Ne, Mg, and Fe to be 
0.1, 1.6, 0.13, 0.22, 0.13 and 0.13 times their solar val- 
ues ([Anders fe Grevessd 119891 ) . respectively, which cor- 
respond to the nebular abundances. The simulated two- 
temperature APEC model spectra were then absorbed 
by an interstellar absorption column and convolved with 
the EPIC-pn response matrices. The resulting spectra 
were then compared to the observed spectrum in the 
0.3-1.3 keV energy range and x 2 statistics are used to 
determine the best-fit models. A minimum of 50 counts 
per bin was required for the spectral fit. The foreground 
absorption (Nn), plasma temperatures (kT±, kT?) with 
l-cr uncertainties, and emission measures (EMi, EM2) 
of the best-fit models are listed in Table [2] The best- 
fit models are overplotted on the background-subtracted 
spectra in FigureO together with the residuals of the fits. 
Multi-temperature models do not provide a substantial 
reduction of the value of the reduced x 2 °f the fit . We 
note that the values of the reduced x 2 differ the most 
from unity for large regions, implying inconsistencies of 
the relative calibrations across the FoV, but the spec- 
tral fits still provide a fair description of the observed 
spectrum. In the following sections we discuss the spec- 
tral fits of the emission from the different morphological 
components of S 308 in more detail. 

Spectral fits using models with varying chemical abun- 
dances of C, N, and O were also attempted, but they 
did not provide any statistical improvement of the fit. In 
particular, models with N/O abundance ratios different 
from those of the nebula resulted in notably w orst quality 
spectral fits. As noted by o ther authors (see IChu et al.l 
120031: iZhekov fe ParkJ 1201 If) , an X-ray-emitting plasma 
with chemical abundances similar to those of the optical 
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Fig. 4. — Composite color picture of the XMM-Newton EPIC image (blue) and CTIO [O III] (green) and Ha (red) images of S 308. The 
apparent X-ray emission outside the optical shell is caused by large fluctuations in the background in regions near the EPIC cameras edge 
where the net exposure is much shorter than at the aimpoint. 



nebulae seems at this moment to be the most adequate 
model for the soft X-ray emission from WR bubbles. 

4.2.1. Properties of the Global X-ray Emission from S 308 

The best-fit to the combined spectrum of the whole 
nebula results in unphysically high values of the hydro- 
gen absorption column density, A^j, well above the range 
[0.2 — 1.05] x 10 21 cm~ 2 implied by the optical extinc- 
tion values derived fr om Balmer decrement of the nebula 
(|Esteban et all [l992l) . The effects of iV H and nitrogen 
abundance on the \ 2 of the spectral fits appear to be 
correlated, i.e.. models with high iVn and low nitrogen 
abundance fit the spectra equally well as models with 



low ./Vh and high nitrogen abundance. If we adopt the 
high absorption column density from the best-fit model 
(JVh >3x 10 21 cm" 2 ) , the e levated nitrogen abundance 
reported bv lChu et al.l (|2003f ) will not be reproduced. As 
the high absorption column density is not supported by 
the optical extinction, in the subsequent spectral fits we 
will adopt a fixed Ah of 6.2 x 10 20 cm -2 that is consis- 
tent with the optical extinction measurements. We note 
that this choice results in an imperfect modeling of the 
spectral features in the 0.3-0.5 keV range, as indicated 
by the S-shaped distribution of residuals in this spectral 
region in Fig. [7] If we allow the value of ATh to float 
during the spectral fit, the improvement of the value of 
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TABLE 2 

Spectral Fits of the Diffuse X-ray Emission of S 308 



Region 


Counts 


ATh 


k 1 \ 


a 


( b 
Ji 




kT2 


a 

EM 2 


f b 

J2 




72/ Jl 




X / ■L'O.r 






[lu cm 


[KCV ] 


[cm J 


r — 2 
[erg cm 


s J 


[KCV J 


[cm J 


[erg cm 


s 1 








S308 


10290±19 


6.2 


0.096±0.002 


7.6x10 s6 


2.7x10" 


12 


1.12±0.22 


1.2X10 55 


2.4x10" 


13 


0.090 


2.01 


(=319.6/159) 


Shell 


7920±23 


6.2 


0.092±0.003 


4.9x10 s6 


1.5x10" 


12 












1.54 


(=233.9/151) 


Interior 


2390±120 


6.2 


0.116±0.011 


8.3x10 s5 


5.6x10" 


13 


0.95 


4.8 x10 s4 


7.5x10" 


14 


0.134 


1.35 


(=188.2/139) 


1NE 


965±10 


6.2 


0.094±0.010 


1.7xl0 58 


5.1x10" 


13 


0.95 


l.lxlO 53 


2.5x10" 


15 


0.005 


0.82 


(=46.9/57) 


1NW 


1000±60 


6.2 


0.102±0.009 


4.7xl0 55 


2.0x10" 


13 


0.95 


4.4X10 50 


1.0x10" 


17 


< 10~ 3 


1.13 


(=57.9/51) 


1SE 


540±50 


6.2 


0.094±0.010 


1.8X10 58 


5.1x10" 


13 


0.95 


l.lxlO 53 


2.5x10" 


15 


0.005 


0.82 


(=46.9/57) 


1SW 


1100±50 


6.2 


0.097±0.012 


l.Ox 10 56 


3.9x 10" 


13 


0.95 


2.3 x 10 54 


5.2x10" 


1 1 


0.135 


1.31 


(=69.3/53) 


2NW 


4620±100 


6.2 


0.095±0.003 


2.2X10 56 


7.4x10" 


13 


0.96±0.21 


4.3 xlO 54 


9.8x10" 


14 


0.130 


1.30 


(=137.9/106) 


3NE 


300±33 


6.2 


0.095±0.023 


3.5X10 55 


1.2x10" 


13 


0.95 


1.4xl0 54 


3.1x10" 


14 


0.262 


0.55 


(=9.4/17) 


3NW 


160±25 


6.2 


0.11±0.04 


5.0x10 s4 


2.5x10" 


14 


0.95 


5.6X10 51 


1.2x10" 


16 


0.005 


0.91 


(=27.4/30) 


5NE 


530±50 


6.2 


0.090±0.015 


7.6X10 55 


2.1x10" 


13 


0.95 


9.0xl0 5:i 


1.8x10" 


14 


0.083 


1.06 


(=44.5/42) 


5SE 


820±60 


6.2 


0.103±0.016 


4.4xl0 55 


1.9x10" 


13 


0.95 


2.2X10 54 


5.1x10" 


14 


0.268 


1.05 


(=54.5/52) 


6NW 


400±50 


6.2 


0.112±0.015 


1.8X10 55 


1.0x10" 


13 












1.42 


(=50.0/35) 


6SW 


210±32 


6.2 


0.12±0.05 


8.3X10 54 


5.7x10" 


14 


0.95 


9.7xl0 53 


1.5x10" 


14 


0.270 


1.69 


( = 27.0/16) 


W99 c 


4560 


35 


0.129 


1.0x10 s6 


6.5x10" 


12 


2.4 


4.3 xlO 55 


1.2x10" 


12 


0.185 


20 (= 


=40/2) 


C03 d 
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0.094±0.009 


8.2x10 s6 


7.2x10" 


12 


u ' ' -0.5 


5.1x10 s4 


1.1x10" 


13 


0.015 


1.02 





"EM = JnldV. 



b Obscrved (absorbed) fluxes for the two-temperature models components in the energy range 0.3-1.3 keV. 

1Wriggd l[T99l. 

IChu et all II2003D . 



> 0.5 



1 Raw spectrum of S 308 
" Background spectrum 



+ + + 



0.5 1 

Energy (keV) 

Fig. 5. — Comparison of the background-unsubtracted raw 
EPIC-pn spectrum of S 308 (black) and scaled EPIC-pn back- 
ground spectrum (red). The Al-K line at ~1.5 keV is an instru- 
mental line. 

the reduced x 2 is negible. 

The parameters of the best-fit model, listed in the 
first line of Table show two plasma components at 
temperatures ^l.lxlO 6 K and ~T.3xl0 7 K with an 
observed flux ratio, /2//1 ^0.09, corresponding to an 
intrinsic flux ratio F2/F1 ^0.06. The total observed flux 
is ~ 3 x 10~ 12 erg cm" 2 s _1 . The intrinsic luminosity at 
a distance of 1.5 kpcQ, after accounting for a fraction of 
~ 1/3 of S 308 which is not included in the source aper- 



1 See discussion about the distance to WR6 of lChu et al.l (2003). 
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Fig. 6. — Background-subtracted blank sky spectra extracted 
from source and background regions equal in detector coordinates 
to the regions 1NW (top) and 1SW (bottom) of S 308. 



tures considered here, is ~ 2xl0 33 erg s _1 . The emission 
measure of the best-fit to the combined spectrum, along 
with the spatial distribution of the X-ray-emitting gas 
in a spherical thick shell with a thickness ~8' and inner 
radius of ~11', implies an average electron density 
n e ~ 0.1 cm" 3 . We note that the quality of the spectral 
fit is not exceptionally good, but more sophisticated fits 
using multi-temperature models failed to improve the 
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E (keV) E (keV) 

Fig. 7. — Background-subtracted XMM-Newton EPIC-pn spectra of S308 corresponding to the 11 individual source regions shown in 
Figure[T] as well as the combined spectra of the entire nebula, its shell or rim (regions #1 and #2) and the interior region of the shell (#3, 
#5, and #6). Each spectrum is overplotted with its best-fit two-temperature APEC model in the energy range 0.3-1.3 keV, assuming fixed 
abundances, while the lower panel displays the residuals of the fit. 



quality of the fit. The proposed 2-T model, providing 
a fair description of the spectral shape, should be 
considered as a first approximate of the hot gas content 
and its physical conditions. 



4.2.2. Northwest Blowout (Region #2) 

The northwest blowout of S308 has the brightest 
X-ray emission, with a surface brightness ~ 2.0 x 
10~ 18 ergcm~ 2 s _1 arcsec~ 2 and its individual spectrum 



has a high signal-to-noise ratio. The spectral shape is 
consistent with those of the shell spectra, with a promi- 
nent 0.43 keV N VI line, a weaker O VII line, and a clear 
detection of X-ray emission to energies of 0.8-1.0 keV. 
The best-fit parameters are rather similar to those of the 
spectrum of the entire nebula, with a marginally lower 
temperature for the hard component. We will adopt this 
value of the hard component temperature for those re- 
gions whose spectra do not have an adequate count num- 
ber to fit this parameter. 
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4.2.3. The Limb-Brightened Shell 

The diffuse X-ray emission from S 308 has a clear limb- 
brightened morphology surrounding a cavity of dimin- 
ished X-ray surface brightness. The emission from re- 
gions at the rim of this shell (1NE, 1NW, 1SE, and 1SW) 
is relatively bright, with an averaged surface brightness 
of the rim ~ 1.2 x 10 -18 ergcm~ 2 s _1 arcsec~ 2 . All indi- 
vidual spectra of the rim regions show the bright 0.43 keV 
N VI emission line and indications of the weaker 0.57 keV 
O VII emission line. The hard component is faint, except 
for the spectrum of region 1SW. The fit to the combined 
spectrum confirms the temperature of the soft compo- 
nent, but it is not possible to provide statistical proof 
of the detection of the hard component. The fits to the 
individual spectra only provide upper limits for this com- 
ponent, except for region 1SW where it seems relatively 
bright. 

4.2.4. The Central Cavity 

The level of X-ray emission from the innermost re- 
gions of the optical shell of S308 is lower than that 
of its edge, with an averaged surface brightness of 5 x 
10~ 19 erg cm -2 s _1 arcsec - , i.e., ~ 2.5-4.0 times fainter 
than the shell and blowout regions. The combined X- 
ray spectrum of the interior regions shown in Figure [7] 
indicates a stronger relative contribution from the hard 
component. This is indeed confirmed by the spectral fit: 
on average, the hard X-ray component has a flux ~ 13% 
that of the soft component. There is a noticeable lack 
of emission from this component in the region 6NW, but 
otherwise the average contributions derived from the in- 
dividual fits are higher than for the spectra of apertures 
on the shell rim. 

4.2.5. Comparison with Previous X-ray Studies 

Table [5] also lists the best-fit parameters of the spectral 
fits to the diff use X -ray emission from S 308 obtained by 
iWriggd (fl999h and IChu et all |200l . It is worthwhile 
discussing some of t he differenc e s wit h these previous 
X-ray analyses. The IChu et all (|2003j ) joint fit of our 
regions 1NW, 2NW, 3NW and 6NW yields very similar 
results to the ones shown in Table [5J For the second 
thermal component, the d erived tem p eratu res from our 
spectral fits and those of IChu et al.l (2003[ ) are consis- 
tent with each other, but IWriggd ( 1996 ) provides a much 
higher temperature for this component. This discrep- 
ancy highlights the difficulty of fitting the hard compo- 
nent using ROSAT PSPC data given its low spectral 
resolution, as well as the very likely contamination of 
the ROSAT PSPC spectrum of S 308 by unresolved hard 
point sources superposed on the diffuse emission. 

5. THE CENTRAL STAR WR6 (HD 50896) 

The central star of S 308 is WR 6 (a.k.a. HD 50896) 
which has a spectral subtype WN4 (jvan der Hucht et al.l 
119881) . The star is detected by the XMM-Newton EPIC 
cameras in the NE and SW pointings of the nebula. 
The EPIC-pn, EPIC-MOS1, and EPIC-MOS2 count 
rates are 160±4 counts ks _1 , 65±2 counts ks _1 , and 
65±2 counts ks _1 , respectively, from the NE observa- 
tion, and 141±4 counts ks _1 , 53±2 counts ks _1 , and 
52±2 counts ks _1 , respectively, from the SW observa- 
tion using a source aperture of 50" in radius. These 



0.0002 



EPIC-pn 



EP1C-MOS1 




EPIC-MOS2 



10 20 30 40 50 

T-T [ks] 

Fig. 8.— EPIC-pn (top), MOS1 (center), and MOS2 (bottom) 
exposure-map-corrected light curves of WR6 in the 0.2—9.0 keV 
energy band. The time is referred to the starting time of the NE 
observation, 2004-03-15T06:45:41 UTC. 



count rates appear to imply that the X-ray luminosity 
of WR6 diminished by 10-20% from the NE to the SW 
observations, which were only ~ 8 hours apar10. We 
note, however, that these count rates are largely affected 
by vignetting due to the offset position of WR 6 on the 
EPIC cameras. Indeed, the light curves shown in Fig- 
ure [51 after accounting for the effects of vignetting, may 
imply the opposite, i.e., that the X-ray flux of WR6 was 
slightly higher in the second (SW) observation than in 
the first (NE) observation. 

In Figure [9] we present the EPIC background- 
subtrac ted spectra fo r the two different epochs. Fol- 
lowing I Skinner et al.l (|2002l ). we have modeled these 
spectra with a two-temperature VAPEC model with 

ini tial abundances set to those shown in Table 1 

of Ivan der Hucht et al.l (|1986D . The fit allowed the 
foreground absorption column density, temperatures, 
and abundances of N, Ne, Mg, Si, and Fe to vary 
(Sk inner et aLll2002f ). Table 3 displays the parameters 
resulting from our best-fit models: absorption column 
densities Nn, plasma temperatures T, normalization pa- 
rameters ^0, observed (absorbed) fluxes /, and intrinsic 
(unabsorbed) fluxes F. Model fits for the spectra from 
the NE and S W observations are lis ted separately along- 
side those from lSkinner et ail (|2002f) for comparison. The 
column density and temperatures of the two components 
are within 1-a of one another among the three different 
mo dels. The observed flu xes are also consistent, although 
the iSkinner et al.l (|2002h flux seems to be a bit lower, 
while our fluxes are closer to the ones derived from the 

2 T he on-axis XMM-Newton images presented bv lSkinner et af] 
pOOlft revealed a faint X -ray source at a radial distance of ~57" 
from WR6. The flux from this source is (2 ± 1)% that of WR6 
and thus it is not likely that fraction of the emission from this 
source entering into the aperture used for WR6 would contribute 
significantly to the observed X-ray variability. 

3 A = 1 X 10~ 14 / nldV / iivd 2 , where d is the distance, n e is the 
electron density, and V the volume in cgs units. 




TABLE 3 
Spectral Fits of HD 50896 



Parameter 










NE Spectrum 


SW Spectrum 


Skinner ct al. (2002) 


N H [xl0 21 cm- 2 ] 












6 4+ 1 - 


5 9+ 1 ' 2 
°- a -0.9 


1 + " 4 
4 - U -0.6 


fcTi [keV] 












n qo+0.03 


0.28 


6+ - 4 
u -°-0.4 


Ai [cm -5 ] 












7.1 xlO" 3 


5.9xl0~ 3 


3.4xl0" 5 


kT 2 [keV] 












2 5+ 1 - 5 


2.48 


°'°-0.6 


A 2 [cm" 5 ] 












1.8xl0~ 4 


2.2xl0" 4 


l.OxlO -5 


X 2 /DoF 












1.10=137.0/124 


1.10=113.2/103 


1.08=234.5/217 


/i (0.2-10 keV) [xlO - 


12 


erg 


cm 


" 2 s" 


1 


1.14 


1.23 


0.97 


/i (2.5-10 keV) [xlO" 


112 


erg 


cm" 


" 2 s- 


\ 


0.33 


0.39 


0.31 


h (0.2-10 keV) [xlO" 


112 


erg cm 


" 2 s- 


L ] 


0.58 


0.70 


0.49 


Fi (0.2-10 keV) [xlO" 


-12 


erf 


; cm" 


" 2 S" 


- 1 ] 


19.8 


17.6 


2.90 


Fj (2.5-10 keV) [xlO" 


-12 


erf 


; cm" 


" 2 S" 


- 1 ] 


0.35 


0.41 


0.33 


F 2 (0.2-10 keV) [xlO" 


-12 


erf 


; cm" 


- 2 s" 


- 1 ] 


1.10 


1.35 


0.78 



Octo ber 1995 ASCA observations of WR 6 (jSkinner et alJ 
119961 ). The total observed fluxes and the observed fluxes 
of the hot thermal component, /2, listed in Table [3] may 
indicate a hardening of the X-ray emission from WR6 
during the last observation. To assess this issue, we per- 
formed statistical evaluation of the lightcurves showed in 
Figure [8] using the HEASOFT task lestats and found no 
significant variations. Thus, WR6 does not show evi- 
dence of variability in time-scales of hours. 

We would like to point out that the absorption 
column density obtained from our best fi ts are in good 
agreem ent with the values obtained from iSkinner et alJ 
(|2002f h which are higher than that used to fit the soft 
X-ray emission from the nebula. Such higher column 
density values are commonly observed towards massive 
stars such as WR stars and are recogniz ed to be caused 



6. DISCUSSION 

The XMM-Newton images and spectra analyzed in the 
previous sections reveal that the hot plasma in S 308 is 
spatially distributed in a thick shell plus the northwest 
blowout, with most emission being attributable to a hot 
plasma at ~l.lx 10 6 K. For an adiabatically shocked stel- 
lar wind, its temperature is determined by the stellar 
wind velocity, kT = 3^,rriB_V^/16, w here n is the mean 
parti cle mass for fully ionized gas (jDyson fc Williams] 
Il997| ). Therefore, the temperature expected for the 

shocked stellar wind of WR 6, with a terminal wind veloc- 
-l 



by absorption at the base o f the wind ( Cassinelli et al 
198H iCorcoran et all 119941 : iNazel 120091 ISkinner et al 



20ld iGosset et alJl201lf T" 



ity of 1700 km s" 1 (|Hamann fc Koesterkdll998f ) and /x > 
1.3 (jvan der Hucht et al.lll986| ). would be T > 8 x 10 7 K, 
in sharp contrast with the observed temperature. The 
same issue has been pointed out for the WR bubble 
NGC6888 by several authors (see lZhekov fc Park)l20iTI . 
and references therein), and it is al so a common issu e in 
planetary nebulae (e.g., NGC6543: IChu et al.ll200ll) . 
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Electron thermal conduction has been proposed as a 
mechanism capable of reducing the temperature of the 
hot plasma in shocked ste llar wind bubbles. Thermal 
conduction was applied by iWeaver et al.1 (|1977l ) to stel- 
lar wind bubbles to produce a self-similar solution for the 
density and temperature structure in bubbles. The soft 
X-ray luminosities predicted by Weaver et al.'s model 
for the Omega Nebula and the Rosette Nebula, accord- 
ing to the stellar wind parameters of their associated 
young clusters (M17 and NGC 2244, respectively), are 
several orders of magnitude higher than t hose observed 
(jTownslev et all [2001 fPunne et all [20031 ). Thus, the 
standard Weaver et al. model for a stellar wind bubble 
with thermal co nduction cann o t be t aken at face value. 
Recent work by iSteffen et al.l ([2008) in the context of 
planetary nebulae, which are produced in a very similar 
manner to WR wind bubbles, has calculated the time- 
dependent radiation-hydrodynamic evolution of plane- 
tary nebula wind bubbles including thermal conduction 
in ID models with spherical symmetry. In these models, 
the cold shell material from the previous AGB superwind 
phase is evaporated into the hot bubble. Saturated con- 
duction was taken into account in these calculations by 
limiting the electron mean free path and it was found 
that thermal conduction was able to lower the tempera- 
ture and raise the density at the edge of the hot bubble 
enough to explain the soft X-ray emission and low X-ray 
luminosities obs erved in some planetary nebulae (e.g., 
IChu et alJ[20fTlh . 

In the case of WR bubbles. iToala fe Arthur! (|2011D pre- 
sented time-dependent 2D radiation-hydrodynamic mod- 
els of the evolution of the CSM around single massive 
stars, including classical and saturated thermal conduc- 
tion. They found that in the absence of a magnetic field, 
thermal conduction does not seem to play an important 
role in shaping WR bubbles, but that models with ther- 
mal conduction have slightly greater soft-band luminosi- 
ties than those without thermal conduction. They sug- 
gested that the morphology of S 308 could result from a 
star with initial mass of 40 M p whose stellar evolutio n 
model includes stellar rotation (jMevnet fc Maedei!l2003f ). 
They obtained that -20,000 yr after the onset of the WR 
phase, the X-ray-emitting gas will present a clump-like 
morphology with an electron density of n e ~ 0.1 cm -3 
inside an optic al (T ~ 10 4 K) shell wit h radius of ~ 9 pc. 

Whereas the lToala fc Arthur! (|2011l ) models reproduce 
the morphology and X-ray luminosity of wind-blown 
bubbles, their simulations predict higher temperatures 
of the hot plasma that result in X-ray spectra that do 
not match the observed spectral shape. This might im- 
ply that additional physical processes must be taken into 
account. An interesting alternative to thermal conduc- 
tion for the apparently low ionization st ate of the plasma 
is non -equilibrium ionization (NEI). iSmith fc Hughes! 
(|2010[ ) calculate the timescales to reach collisional ion- 
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ization equilibrium (CIE) for ionized plasmas and their 
results suggest that, for values of the parameters rele- 
vant to S 3 08 (derived electron d ensity and time in the 
WR stage: iToala fc Arthurll20ll . the CIE assumptions 
may not hold. These ideas will be pursued in subsequent 
works. 

7. SUMMARY AND CONCLUSIONS 

We present XMM-Newton observations of three fields 
of the WR bubble S 308 which, in conjunctio n with the 
observ ation of its NW quadrant presented by IChu et ah! 
(|2003h , map most of the nebula except for its southern- 
most section. We have used these observations to study 
the spatial distribution of the X-ray-emitting material 
within this bubble, to derive global values for its phys- 
ical conditions (T e , n e ), and to search for their spatial 
variations among different morphological components of 
the nebula. 

The X-ray emission from S 308 is found to have a 
limb-brightened morphology, with a shell thickness 5'- 
8', and extend to the northwest blowout region. The 
X-ray-emitting shell is notably larger along the SE-NW 
direction than along the SW-NE direction, and it is al- 
ways confined by the optical shell of ionized material. 
The X-ray surface brightness decreases notably from the 
blowout region and the western rim shell to the shell inte- 
rior, where the X-ray emission falls to background levels. 
The western quadrants are also brighter than the eastern 
quadrants. 

The X-ray emission from S 308 shows prominent emis- 
sion from the He-like triplet of N VI at 0.43 keV and 
fainter emission of the O VII 0.57 keV triplet, and de- 
clines towards high energies, with a faint tail up to 1 keV. 
This spectrum can be described by a two-temperature 
optically thin plasma emission model with temperatures 
~ 1.1 x 10 6 K and ~ 1.3 x 10 7 K. The latter component 
is notably fainter than the former by at least a factor 
of — 6. There is an appreciable difference in the rela- 
tive contributions of the hot component to the X-ray- 
emitting gas between the rim and the nebula interior, of 
which the latter has a higher contribution from the hard 
component. The total X-ray luminosity is estimated to 
be — 2 x 10 33 erg s _1 for a distance of 1.5 kpc. 
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